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ABSTRACT 

We differentiate between tfie metal enricfiment of tlie gas in virialized mini- 
halos and that of the intergalactic medium at high redshift, pertinent to cos- 
mological reionization, with the initial expectation that gas in the high density 
regions within formed dark matter halos may be more robust thus resistant to 
mixing with lower density intergalactic medium. Using detailed hydrodynamic 
simulations of gas clouds in minihalos subject to destructive processes associated 
with the encompassing intergalactic shocks carrying metal-enriched gas, we find, 
as an example, that, for realistic shocks of velocities of 10 — lOOkm/s, more than 
(90%, 65%) of the high density gas with p > 500p;, inside a minihalo virialized 
at z = 10 of mass (10'', 10^)Mq remains at a metallicity lower than 3% of that 
of the intergalactic medium by redshift z = 6. It may be expected that the 
high density gas in minihalos becomes fuel for subsequent star formation, when 
they are incorporated into larger halos where efficient atomic cooling can induce 
gas condensation hence star formation. Since minihalos virialize at high redshift 
when the universe is not expected to have been significantly reionized, the im- 
plication is that gas in virialized minihalos may provide an abundant reservoir 
of primordial gas to possibly allow for the formation of Population-Ill metal-free 
stars to extend to much lower redshift than otherwise expected based on the 
enrichment of intergalactic medium. 

Subject headings: supernovae: general — galaxies: formation — intergalactic 
medium — cosmology: theory — early universe 



Introduction 



Recent observa tions of the high redshift (z > 6) quasar spectra frora the Sloan Digital 



Sky Survey (SDSS; iFan et al. II2OOII : Becker et al.ll200ll : iBarkana I I2OO2I : ICen fc McDonald 
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20021 ) and the cosmic micro wave background flu ctuations from the Wilkinson Microwave 
Anisotropy Probe (WMAP; ISpergel et aL 1 12006| ) combine to paint a comphcated and yet 
uncertain reionization picture. What is strongly suggested, however, is that the reionization 
process began at an early time, probably at z > 15. If stars are primarily responsible for pro- 
ducing the ionizing photons, the overall reionization picture would depend significantly on 



how and when the transition from me tal-free Popula t ion III (Pop II 



tion 1 1 (Pop II) stars occurs (see. e.g..lOmukai 



2OO2I : iMackev et al. II2003I : ISchneider et al. 
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.2001al : ISchneider et al. 
2003h. especiallv if Pop 



III starformation history is more extended than usually thought (IWyithe &: Cen 1 120071 ). 



This is based on the finding t 


lat Pop '. 


II stars may be much more massive and hotter than 


Pod II stars ( Carr et al. 1 1984 


Larson 


Il998: Abel et al. 2000: Hernandez & Ferrara 2001: 


Bromm et al. 2001a 


b: Nakamura & Umemura II2OOI: Bromm et al. 2002: Omukai & Palla 
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2003) and hence efficient producers of ionizing photons of hvdrogen and 



helium. 

The metallicity of the star-forming gas plays several important roles in the physics of 
first stars. First, the transition from Pop III to Pop II is facilitated by the pre sence of a 
small amount of metals, in particular, oxygen and carbon (IBromm fc Loeb 1120031 ). Thus, it 
is the amo unt of C and O, not necessarily the total amount of "metals" , that determines the 
transition (IFang fc Cen 112004 ). The yield patterns for (non-rotating) stars with mass in the 
range of 140 — 260Mq that explode via the pair-instability supernovae (PISN) and regular 
type II SNe are different. In PISN case the supernova ejecta is enriched by a-elements, 
whereas the major products of SNII are hydrogen and helium with a sm all amount of heavy 



elements (see, e.g.. IWoosley fc Weaver Ill995l : iHeger fc Woosley II2OO2I ). Consequently, the 



trans ition from Pop III t o Pop II stars may occur at different times, depending on the IMF 
(e.g., iFang fc Cen"fl2004h . 



Second, while the ionizing photon production efficiency depends only relatively weakly 



on th e exact IMF, as long as the stars are more massive than ~ IOMq (e.g.. lTumlinson et al. 
2OO4J ). its dependence on metallicity is strong, because the effective temperature of the stellar 
photosphere depends sensitively on the opacity hence metallicity of the stellar atmosphere. 
The amount of metals produced depends on the IMF. For example, in the most extreme case 
where all Pop-Ill stars are more massive than, say, > 27OM0, these stars may conclude by 
implosions to intermediate-mass black holes without giving out much metals to the surround- 
ings. However, exactly how massive Pop HI stars are is uncertain. While simulations have 
suggested that Pop HI stars may be mor e massive than 100 Mf7> ( "very massive star", VMS; 
Abel et al. I I2OO0I : iBromm et allboOlah . iTan fc McKeel JiooJ find that stellar feedback 
processes may limit the mass of th e Pop HI s t ars to the r ange 30 — 100 M(?^. Observat ionally, 
the VMS picture is advocated by lOh et al. I (1200 1[ ) and iQian &: Wasserburg I (120021 ). based 
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on an analysis of metal yield patterns fro m pair-instability superno v a (PISN) explosion o f 



VMS progenitors (iHeger fc Wooslev II20021). iTumlinson et a. 



f|2004[ ). lDaigne et al. I (120041 ) 



Umeda & Nomoto (2003,2005) and IVenkatesan &: Truran I (120031 ) . on the other hand, argue 
that the general pattern in metal-poor halo stars, in the Lya forest and cosmic star forma- 
tion history, is more consistent with the yield pattern of Type II supernovae (SNII) perhaps 
with a lower cutoff of lOM©. 

Clearly, the metallicity of gas out of which stars are formed is critically important. The 
conventional picture that is often adopted goes as follows: formed stars will eject metals 
into the IGM and eventually raise the metallicity of the IGM to above the threshold for the 
Pop-Ill to Pop-II transition. A somewhat refined version of this takes into account that the 
meta llicity enrichment process of the IGM is unlikely to be synchronous for different regions 
(e.g., Furlanetto &: Loebl 2005 ). Here, we point out a possibly large difference between the 



metallicity of the IGM an d the metallicity of th e gas in minihalos. Since minihalos collapse 
at very high redshift (e.g.. IWyithe fc Gen 1120071 ). the large amount of dense gas in minihalos 
thus may provide a primary fuel for subsequent star formation, when eventually they are 
incorporated into large systems where efficient atomic cooling allows gas to condense to form 
stars. 

To quantify this possible difference between the metallicity of minihalo gas and that of 
the IGM, we study the stability and metal enrichment of minihalos subject to metal-rich 
Shockwaves launched by supernovae explosions from large galaxies. We will treat an ideal- 
ized situation where a minihalo is subject to shock waves enriched with a chosen metallicity, 
and we investigate how gas inside it may be contaminated by metals. We will assume that 
there has been no star formation hence no self-metal-enrichment in minihalos, because of 
the lack of adequate coolants; molecular hydrogen is assumed to have long been destroyed 
by Lyman- Werner photons produced by earlier stars elsewhere. Because gas in minihalos 
is significantly overdense compared to the IGM and is bounded by the gravitational poten- 
tial wells produced by the dark matter halos, mixing of metals into the gas in minihalos 
by metal-rich outfiows from star-forming galaxies should be expected to be different from 
that of the IGM. As we will show, the process of mi xing of metal-rich ou tflows with the 



gas in minihalos is quite incomp l ete. Several authors (IMurray et al. I Il993l . M93 hereafter; 



Klein et al. 1 1 1994 iDinge 1 119971 : iMiniati et al. 1 119971 ) have addressed the problem of the 



stability of a non-self-gravitating gas cloud moving at the sound speed of the background 
medium, which is equivalent to a shockwave sweeping the gas cloud. They have found that 
the cloud gets disrupted after a time comparable to the dynamical time of the cloud. Here, 
we are interested in the self-gravitating case. In particular, we are interested in minihalos 
that are gravitationally dominated by their dark matter content, and with no cooling. A 
very similar case was already studied by M93, in the context of a two-phase medium, using 
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2-D simulations. In this work we employ 3-d hydrodynamical simulations to study this prob- 
lem. We simulated halos of mass IO^Mq and IO^Mq subject to Shockwaves with velocities 
of 10, 30, 100, and 300 km/s. For the slowest cases of 10 and 30 km/s the halos are quite 
stable and the gas inside the virial radius of the halos remains fairly uncontaminated after 
many dynamical times. Only for shock velocities of 100 and 300 km/s the halos start to be 
unstable, loosing significant fraction of their gas, and getting substantially enriched in their 
inner regions. 

The paper is organized a follows. In §2 we specify the physical model for the minihalos 
and Shockwaves, and describe some technical specification for the code we use. §3 presents 
our results, followed by conclusions in §4. 



2. Description of the model 



We analyze the metal enrichment of gas in spherical minihalos with total virial masses 
of 10^ and lO^M©, whose virial temperatures are 710 K and 3295 K, respectively, at z = 10. 
Initially, the gas in minihalos is assumed to have zero metallicity. Then, the minihalo is 
exposed to an IGM sweeping through at a velocity of K and metallicity Zjgm, and we 
quantify the evolution of the metallicity of the gas inside the minihalo. We study four cases 
with Vg = 10, 30, 100, and 300 km/s for each of the two choices of minihalo masses. 

The gravitational potential of a halo is determined by their dark matter and assumed 
not to change. Th e density of a yirializ ed dark matter halo as a function of radius, r, is 
given by the NFW (INavarro et al. 1119971 ) density profile: 



PDAiir) 



Pcrit^c 



where pcrit = 3H{zY /SnG is the critical density of the universe at redshift z, Sc = 200c^/3 m(c), 
and u = r/rs- The characteristic radius is defined in terms of the concentration parameter 
of the halo, c, that is a function of the halo mass and the redshift, and the virial radius, r^r- 
The virial radius is defined in terms of the halo mass, M^, by (47r/3)r^j^200pcrjt = M^, and 
the function m{u) = ln(l + u) — u/(l + u). Fo r the concentration parameter we adopt the 
fitting formula provided by iDolag et al. I (12004J ) : 



9.59 



M, 



H 



-0.102 



1 + z'lO^^h-^Mr. 







based on computations of a ACDM cosmological model with 
and (jg = 0.9. 



(2) 



0.3, = 0.7, = 0.045, 
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Since the gravitational potential, 0, is determined by the dark matter content of the 
minihalos, it will be given by: 



>ln(i+u) 



-A-nG ScPcrit f1 — - — — foTU<d 

crcrti, s __ u _ — ^g'j 

for u > d. 



where Vgd is the radius at which the dark matter density of the halo equals the mean density 
of the universe. 

For the gas we used a. X = 0.76, Y = 0.24, and Z = composition. Its density and 
temperature profiles will be determined by assuming Pg = pgksTg/ {nipfi) a pj, where 7 is 
the polytropic index. Then we can write 

Pgir) = pcViu) and Tg{r) = Tcyiuf^^ , (4) 

where pc and Tc are the density and temperature at the center of the halo. Since the halos 
are in hydrodynamic equilibrium we find, 

( ^ 7-l GAfHmp^nn(l+n) _ ^\ for M < 

We have three free parameters in our gas profile: pc? T'^, and 7. The central density will be 
chosen such that at r = Vmr the ratio between dark and baryonic matter densities is equal 
to fim/f^b. The values for Tc and 7 are determined considering that, according to many 
hydrodynamic simulations, in the outer part of the h alos the gas density profile traces the 
dark matter density profile (IKomatsu fc Seljak Il200ll ). i.e., 

d\Ty{pDM) ^ 

dlnipg) ■ ^ ' 



So Tr and 7 were determined using the fitting formulae provided by iKomatsu fc Seljak 
(I2OOII ) that satisfy equation within a range c/2 < u < 2c: 

Tc = ^^gg^(0.00676(c- 6.5)2 + 0.206(c- 6.5) + 2.48) 

7 = 1.15 + 0.01(c- 6.5). ^ '* 

As mentioned above, the minihalos will be shocked by a continuous, hot, metal-rich shock- 
wave at a temperature Tg = ^j^^, a mean density ps = 4:pi„ and a metallicity Zjgm- We 
also introduce fluctuations in space and time in the density of the shockwave, i.e., 

2txv 271 z 27ct 

= 4pb(l + Asin(^ + 03,)sin(— + 0,)sin(— + 0^)), (8) 
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where y and z represent the two spatial coordinates perpendicular to x, the direction of 
propagation of the shockwave, t is the time, A and A represent the amplitude and the 
length of the fluctuation, and (pi corresponds to an arbitrary phase in the coordinate i. 
The spatial phases were randomly chosen every time that t/T became an integer number, 
where T = X/Vg. What will be appropriate values for A and A? At the redshifts of interest 
(z ~ 6 — 10), large atomic-cooling halos start to become nonlinear, meaning that the density 
variance is of order unity on the mass scales of 10^ — IO^Mq, corresponding length scales of 
~ 0.1 Mpc. By definition, when a certain mass scale M become nonlinear, A{M) ~ 1. We 
have experimented with values oi A = 0.3 — 0.9 and A = 0.003, 0.01, and 0.03 Mpc/h in 
comoving units. Our results turned out to be nearly independent of the values for A and A 
in the ranges of relevance. 

Each simulation starts dX z = 10, when the IGM shockwave enters the left face of our 
simulation cube. We do not attempt to vary the background density with time, aside from 
the variation imposed (see Equation ([8]) above). We expect that, if the background density 
were allowed to decrease with time, the metal enrichment of the gas in minihalos may be 
reduced. When there is a need to indicate a redshift during the evolutionary phase of a 
minihalo, we translate the elapse d time since ^ = 10 t o a certain redshift, using the standard 



cosmo logical model parameters (ISpergel et al. II2006I ) 



We use the TVD hydrodynamics code (ICen et al. II2003I ) to perform the described sim- 



ulations. The size of the boxes is chosen such that at the border of the box the gas density 
of the halo is equal to the mean baryonic density of the universe. So the comoving size of 
the boxes is 0.0191 and 0.0457 Mpc/h for Mh = 10^ and lO^M©, respectively. For most of 
the simulations we use 256'^ cells for each simulation. Our results seem to be convergent to 
a few percent accuracy as will be shown at the end of the next section. 



3. The results 

In this section we analyze our results for the stability and chemical evolution from z = 10 
to 2; = 6 for the two halos considered. Although for parameters of the range considered we 
observe different levels of instability and mixing, in all cases the gas in the inner region of the 
halos remains substantially less metallic than the IGM. For Mh = 10'' Mq and Vs = 10 km/s 
almost all the mass at density higher than the virial density, p^r (= 49.3 pi, for = lO^M©) 
ai z = 6 has Z < O.OSZjgm, whereas for Mh = 10^ Mq and Vg = 300 km/s most of the mass 
at p > pmr (= 40.3 pb for Mh = IO^Mq) remains with Z < O.SZjgm- 

Figures ([T]) and ([2]) show density, metallicity and velocity of the gas on a shce through the 
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Fig. 1. — The density and metallicity in a slice through the center of the box for Vg = 30 km/s 
and Mh = 10'' M©. The upper and bottom plots show density and metallicity, respectively. 
The left plots correspond to 2; = 9 and the right ones to 2; = 6. The velocity field along with 
contours of p = pmr (solid line) and p = 500pb (dotted line) are depicted in all the plots. 
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Fig. 2. — The density and metallicity on a slice through the center of the box for Vg = 
100 km/s and Mh = IO^Mq. The upper and bottom plots show density and metallicity, 
respectively. The left plots correspond to 2; = 9 and the right ones to 2; = 6. The velocity 
field along with contours of p = p^r (solid line) and p = 500pb (dotted line) are depicted in 
all the plots. 



- 9 - 



center of the halo that is perpendicular to the shockwave front, at z = 9 and z = 6. Perhaps 
the most noticeable is that the gas cloud inside the minihalo is able to withstand significant 
Shockwaves and to reside inside the halo gravitational potential well for an extended period 
of time. The gravitational potential well of the dark matter halo is more "steady" than a 
pure self-gravitating gas cloud, since the dark matter is unaffected by gasdynamic processes 
at the zero-th order, in agreement with M93. 

However, the mixing between the primordial minihalo gas and metal-enriched IGM due 
to hydrodynamical instabilities is apparent. First, as seen in the left panels in Figure ([2]), the 
Richtmyer-Meshkov instability seems most apparent when the interface between the sweeping 
IGM and the minihalo gas cloud is first being accelerated by the shock moving from left to 
right. Subsequently, with the build-up of a smoother and larger density transition region 
on the left side of the halo gas cloud and reduced shock strengths, the Richtmyer-Meshkov 
instability progressively abates. Second, as can be seen in the right panels in Figures ([1]) 
and ([2]), the Kelvin- Helmholtz instability provides an efficient mechanism to mix gas in the 
shearing regions at the outer part of the minihalos. The fact that the density peak largely 
remains at the center of the dark matter halo over the extended period of time while the outer 
layers become mixed with the IGM suggests that mixing due to hydrodynamic instabilities 
plays the dominant role, whereas ram-pressure stripping is sub-dominant, at least for these 
two cases and during the displayed time interval. Nevertheless, the central regions of the 
minihalo gas clouds are significantly contaminated with metals at later times (right panels in 
Figures [U [2]) . We will later show in Figure ([7]) our convergence test of the results, suggesting 
that our numerical resolution appears to be adequate to properly simulate hydrodynamic 
instabilities involved. 

We will now turn to more quantitative results, focusing on the metal enrichment of gas 
inside minihalos by the IGM shocks. Figures ([MH]) show the evolution of the amount of mass 
at p > 500pfc and p > p^r respectively, in units of its corresponding value at the beginning 
of the simulation at z = 10, that is metal-enriched to various levels with Z < oZjgm with 
a = (1,0.3,0.1,0.03). Figures §Ml show cases with Vs = (10, 30, 100, 300)km/s. 

From Figure with Vg = lOkm/s we see that for = IO^Mq only 5% of the gas 
is contaminated to Z > O.OSZjgm ^ot p > p^r, and for a IO^Mq halo there is practically 
no gas with Z > O.OSZjgm in that density range after about 11 dynamic time (by z ~ 6). 
For p > SOOpfe, there is no gas with Z larger than 0.03Z/g'm even for the Mh = lO^M© 
halo. It is interesting to note that for this velocity the amount of gas at the two ranges of 
density considered actually increases instead of decreasing. This is due to the compression 
produced on the halo by the gas from the shock. We also observe the acoustic oscillations 
in the amount of gas due to this compression. (The fact that acoustic oscillations for the 
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Fig. 3. — The evolution of the mass of gas in different ranges of density and metalhcity 
for a velocity shock of 10 km/sec. The masses considered are 10^ and IO^Mq, represented 
by thin and thick lines, respectively. The left plot takes density range p > 50Qph, where 
Pb is the mean baryon density of the universe at z—10. The right plot takes density range 
P > Pvir, where pmr is the density at the virial radius of the halo (~ 40.3pb and 49.3pb for 
a 10^ and IO^Mq halos, respectively). The metalhcity ranges considered are Z / Zjq]^^ < o;, 
where a — 1,0.3,0.1,0.03 are represented by solid, dotted, dashed, and dot-dashed lines, 
respectively. The dynamical time is tdyn — {SOOirGflrnPcrit)'^ ■ 
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Fig. 4. — The evolution of the mass of gas in different ranges of density and metalhcity 
for a velocity shock of 30 km/sec. The masses considered are 10^ and IO^Mq, represented 
by thin and thick lines, respectively. The left plot takes density range p > 5QQpb, where 
Pb is the mean baryon density of the universe at z—10. The right plot takes density range 
P > Pvir, where pmr is the density at the virial radius of the halo (~ 40.3pb and 49.3pb for 
a 10^ and IO^Mq halos, respectively). The metalhcity ranges considered are Z / Zjq]^^ < o;, 
where a — 1,0.3,0.1,0.03 are represented by solid, dotted, dashed, and dot-dashed lines, 
respectively. The dynamical time is tdyn — {SOOirGflrnPcrit)'^ ■ 
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Fig. 5. — The evolution of the mass of gas in different ranges of density and metalhcity 
for a velocity shock of 100 km/sec. The masses considered are 10^ and lO^M©, represented 
by thin and thick lines, respectively. The left plot takes density range p > 500pft, where 
Pb is the mean baryon density of the universe at z=10. The right plot takes density range 
P > Pvir, where pmr is the density at the virial radius of the halo (~ 40.3p6 and 49.3p6 for 
a 10^ and IO^Mq halos, respectively). The metalhcity ranges considered are Z / Zjq]^^ < en, 
where a — 1,0.3,0.1,0.03 are represented by solid, dotted, dashed, and dot-dashed lines, 
respectively. The dynamical time is tdyn — {SOOirGflrnPcrit)'^ ■ 
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Fig. 6. — The evolution of the mass of gas in different ranges of density and metalhcity 
for a velocity shock of 300 km/sec. The masses considered are 10^ and lO^M©, represented 
by thin and thick lines, respectively. The left plot takes density range p > 500pft, where 
Pb is the mean baryon density of the universe at z—10. The right plot takes density range 
P > Pvir, where pmr is the density at the virial radius of the halo (~ 40.3p6 and 49.3p6 for 
a 10^ and IO^Mq halos, respectively). The metalhcity ranges considered are Z / Zjq]^^ < en, 
where a — 1,0.3,0.1,0.03 are represented by solid, dotted, dashed, and dot-dashed lines, 
respectively. The dynamical time is tdyn — {SOOirGflrnPcrit)'^ ■ 
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10^ Mq halo start earlier is just due to the smaller simulation box size.) Figure (jlj) shows 
the case Vs = 30 km/s. Here we see that for p > SOOpt again there is not gas that gets more 
metal-rich than Z = O.OSZjgm by z ~ 6. For p > p^r, only ~ 5% of the gas ends up with 
Z > O.OSZjGM for a 10''' Mq halo. However, for = IO^Mq, ~ 5% of the gas mass reaches 

Z > O.lZiGM- 

For Vs = 100 and 300 km/s (Figures [5], [6]) the stripping of the outer parts of the halo 
becomes more important and we start to see that the amount of metal-free gas for the density 
ranges considered starts to decrease significantly for two reasons. First, the halo is losing 
a significant amount of its mass and, therefore, its global structure is being modified. So 
we observe a decrease in the total amount of mass for p > 500pb and p > pmr. Second, 
this stripping put into contact the IGM gas with the innermost part of the halo, moving 
the mixing layer inward and increasing the efficiency of the mixing to higher density regions 
in the halo. For Vg = 100 km/s, at p > SOOpt there is not significant mixing but just a 
small overall reduction of the mass. On the other hand, for p > p^r the total decrease of 
mass starts to be significant reaching even ~ 50% for Mh = IO^Mq, and the amount of 
gas purer than 0.03Z/gm is only ~ 30% and ~ 50% of the original counterparts at z = 10 
for Mh = IO^Mq and Mh = lO^M©, respectively. For K = 300 km/s, at p > 500p6 we 
observe significant reduction of the overall mas, especially for Mh = IO^Mq where the mass 
is reduced to ~ 30% percent of its original value. We can see that the mixing itself does not 
play a very significant role at these densities, with practically no difference between the total 
mass and the mass of gas with Z < 0.03Zjgm for Mh = IO'^Mq. The same thing happens 
with Mh = IO^Mq, but in this case for Z < O.IZjgm- For p > p^ir we see that, along with 
the total reduction of mass, there is substantial enrichment of the gas. The mass of gas with 
Z < 0.03Z/GAf is only ~ 3% and ~ 25% of its value at z = 10 for Mh = 10^ and lO^M©, 
respectively. 

To provide a convenient numerical form we summarize in Tables ([T]) and ([2]) the masses of 
gas in the different ranges of density and metallicity for the halo masses and shock velocities 
considered at the relevant redshifts z = 7 and z = 6, respectively. 

The sensitivities of the gas cloud disruption on shock velocity may be understood in the 
context of instability analysis by M93. M93 show that, when the parameter rj, defined as 

gDRci , . 

is above unity, the cloud is stable up to many dynamical times, where D is the density ratio 
of the gas cloud to the background gas, Rd is the radius of the gas cloud and g is the surface 
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gravity. Numerically, 



r^(r) = ( )2/3.1±fwi^N-4.7 /lo^ 

'^^ ^22km/s' ^WmJ ^ 11 '^Mh' ^ ' 



wher e we have assumed that the density slope near the virial radius is —2.4 (jNavarro et al. 



19971 ): Mr is the mass with radius r, z is redshift. Equation (10) suggests that for Vs < 
25km/s, the gas cloud in minihalo of mass Mh = lO^M© at z ~ 10 is generally quite stable, 
in agreement with our results. For Vs = 300km/s and Mh = lO^M© one obtains rj = 0.005 
at z = 10 and Mr = Mh, suggesting that the outskirts of the minihalo gas cloud would be 
disrupted on the order of a dynamic time, consistent with our results. For Vg = lOOkm/s and 
Mh = IO^Mq we find r] = 0.23 at z = 10 and M, = Mh; M93 find that at r] = 0.25, the gas 
mass loss is still relatively small over many dynamic times, consistent with our simulations. 



Mori et al. I (120021 ) show, in simulations of propagation of supernova blastwaves from 
10^/i~^Mq galaxies at z = 9, that after more than a hundred million years the relative filling 
factor for regions being swept by shocks of velocities larger than U = (10, 30, 100)km/s is 
roughly (100%, 35%, 10%). We expect the velocities to be still smaller at the higher redshifts 
of concern here, due to enhanced cooling and larger Hubble velocity. Therefore, in real 
cosmological setting, combined with our findings, we expect that a large fraction of the gas 
already virialized with minihalos will be largely unaffected by metal- carrying blastwaves and 
remain metal-free to modern redshift, possibly as low as z = 5 — 6, when gas in minihalos 
may be photo-evaporated globally. 

It is prudent to check the convergence of computed results. We have performed addi- 
tional simulations with 128^ and 512'^ grid points. We show in Figure ([7]) an example of these 
convergence tests we have done. We see that, while the difference between the 128^ and 256'^ 
cases can amount up to tens of percent at late times (say, t/tdyn > 5), the difference between 
the 256^ and 512'^ cases is dramatically reduced and is at a level of a few percent even at 
very late times {t/tdyn > 10). It is instructive to notice that the tendency is to decrease 
the level of mixing as we increase the resolution. So our results must be interpreted as an 
upper limit in the metal enrichment of minihalos by Shockwaves, with an accuracy of a few 
percent. 



4. Conclusions 

It is frequently assumed that the metallicity of the intergalactic medium is the primary 
determinant of the epoch of the transition from Pop-Ill to Pop-II stars. We wish to point 
out a potentially large difference between the metallicity of the intergalactic medium and 
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Fig. 7. — Simulations with different resolutions were run to test convergence. Results of the 
evolution of the mass of gas in different ranges of density and metallicity using 128^, 256^, 
and 512^ cells are presented in thick, intermediate, and thin lines, respectively. In the three 
cases a 10'^ Mq halo with a shock velocity of 100 km/sec was simulated. The left plot takes 
density range p > 500pf,, where ph is the mean baryon density of the universe at z=10. The 
right plot takes density range p > pmr, where p^r ~ 49.3p5 is the density at the virial radius 
of the halo. The metallicities considered are Z/Zjgm < oi, where a = 1,0.3,0.1,0.03 are 
represented by solid, dotted, dashed, and dot-dashed lines, respectively. The dynamical time 

is tdyn = {SOOTrGQrnPcrity^- 
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Table 1: Mass of gas in different ranges of density and metallicity at t = 7.07 tdyn (z = 
7), where tdyn — (8007rGQ^Pcrit)~^. Results are presented for the two kind of halos (10^ 
and IO^Mq) and the four shock velocities {Vg — 10, 30, 100, and 300 km/sec) considered. 
The density ranges considered are p > 500pb, where pb is the mean baryon density of the 
universe at z=10 and p > p^r, where pyir is the density at the virial radius of the halo 
40.3pb and A9.3pb for the 10^ and IO^Mq halos, respectively). The mass of gas is measured 
in units of the total mass at a given range of density when t = (z=10). The metallicity 
ranges considered are Z/Zigm < 0.3, 0.1, and 0.03. 



z^7 

t — 7.07 tdyn 


Mh = IO^Mq 


Mh = WMq 


P > P500 


P Pvir 


P > PbOO 


P Pvir 


Vg — lOkm/sec 


Z/ZiGM < 0.3 


1.069 


1.015 


1.294 


1.184 


Z/ZiGM < 0.1 


1.069 


1.013 


1.294 


1.184 


Z/ZrGM < 0.03 


1.069 


1.003 


1.294 


1.184 


Vs = 30km/sec 


Z/ZiGM < 0.3 


0.897 


1.096 


1.053 


0.977 


Z/ZiGM < 0.1 


0.897 


1.072 


1.053 


0.972 


Z/ZiGM < 0.03 


0.897 


0.936 


1.053 


0.953 


Vs = lOOkm/sec 


Z/ZrGM < 0.3 


0.846 


0.947 


0.997 


1.111 


Z/ZrGM< 0.1 


0.846 


0.839 


0.997 


1.055 


Z/ZiGM < 0.03 


0.842 


0.630 


0.997 


0.868 


Vs = 300km/sec 


Z/ZiGM < 0.3 


0.574 


0.765 


0.785 


1.003 


Z/ZiGM < 0.1 


0.573 


0.505 


0.785 


0.823 


Z/Zja_u < 0.03 


0.150 


0.235 


0.765 


0.557 
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Table 2: Mass of gas in different ranges of density and metallicity at t = 11.2 tdyn (z = 
6), where tdyn — (8007rGQ^Pcrit)~^. Results are presented for the two kind of halos (10^ 
and IO^Mq) and the four shock velocities {Vg — 10, 30, 100, and 300 km/sec) considered. 
The density ranges considered are p > 500pb, where pb is the mean baryon density of the 
universe at z=10 and p > p^r, where pyir is the density at the virial radius of the halo 
40.3pb and 49. Sp^ for the 10^ and IO^Mq halos, respectively). The mass of gas is measured 
in units of the total mass at a given range of density when t = (z=10). The metallicity 
ranges considered are Z/Zigm < 0.3, 0.1, and 0.03. 



z = 6 

t — 11.2 tdyn 




Mh = WMq 


P > P500 


P Pvir 


P > PbOO 


P ^ Pvir 


Vg — lOkm/sec 


Z/ZiGM < 0.3 


1.095 


1.041 


1.294 


1.155 


Z/ZiGM < 0.1 


1.095 


1.024 


1.294 


1.155 


Z/ZrGM < 0.03 


1.095 


1.009 


1.294 


1.155 


Vs = 30km/sec 


Z/ZiGM < 0.3 


0.926 


0.944 


1.054 


0.990 


Z/ZiGM < 0.1 


0.926 


0.905 


1.054 


0.967 


Z/ZiGM < 0.03 


0.926 


0.754 


1.054 


0.948 


Vs = lOOkm/sec 


Z/ZrGM < 0.3 


0.680 


0.556 


0.906 


0.811 


Z/ZrGM< 0.1 


0.680 


0.464 


0.906 


0.700 


Z/ZiGM < 0.03 


0.667 


0.347 


0.906 


0.508 


Vs = 300km/sec 


Z/ZiGM < 0.3 


0.288 


0.313 


0.622 


0.553 


Z/ZiGM < 0.1 


0.286 


0.160 


0.622 


0.381 


Z/Zja_u < 0.03 


0.087 


0.031 


0.587 


0.211 
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the metallicity of the gas in minihalos. Utihzing hydro dynamic simulations of gas clouds 
in minihalos subject to destructive processes associated with the encompassing intergalactic 
shocks carrying metal-enriched gas, we find that a large fraction of gas in virialized minihalos 
remains at a metallicity much lower than that of the intergalactic medium. For example, 
for realistic shocks of velocities of 10 — lOOkm/s, more than (90%, 65%) of the high density 
gas with p > 500pb inside a minihalo virialized at 2; = 10 of mass (10^, 1O^)M0 remains at 
a metallicity lower than 3% of that of the intergalactic medium by redshift z = 6, under 
the harsh condition that the minihalo is exposed to Shockwaves continuously from 2; = 10 to 
z = 6. 

In the standard cosmological model, if large halos with efficient atomic cooling are 
responsible for producing most of the reionizing photons, smaller minihalos virialize before 
the universe is significantly reionized. Thus, gas in virialized minihalos may provide an 
abundant reservoir of primordial gas to possibly allow for the formation of Population-Ill 
metal-free stars to extend to much lower redshift than expected otherwise based on the 
enrichment of intergalactic medium. 

A related issue that is not addressed here concerns the fate of the gas inside minihalos 



of the photo-evapo ration of gas in mini halos (IBarkana fc Loeb 



Ciardi et al. I l2006l ) may be ~ lOOMyrs (IShapiro et al. II2004J ): the timescale may be still 



2002 




Iliev et al. 


2005 



lo nger at higher r edshif ts [z > 10) and/or at lower ionizing fiuxes than used in the work 



of IShapiro et al. I (120041 ). It may be that a full understanding requires detailed calculations 



that incorporate both radiative transfer and metal-enrichment processes. 
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